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ABSTRACT 



m 

>- ■ It has been suggested that quasi-periodic oscillations of accreting X-ray 

in 
in 

' axisymmetric linear perturbations to an isentropic, isothermal, unmagnetized 



O 



sources may relate to the modes named in the title. We consider non- 



thin accretion disk. The radial wave equation, in which the number of vertical 



^ ■ nodes (n) appears as a separation constant, admits a wave-action current that 

■ is conserved except, in some cases, at corotation. Waves without vertical nodes 
Q^' amplify when reflected by a barrier near corotation. Their action is conserved. 
Q ■ As was previously known, this amplification allows the n = modes to be unsta- 
^ • ble under appropriate boundary conditions. In contrast, we find that waves with 

■ n > are strongly absorbed at corotation rather than amplified; their action is 
not conserved. Therefore, non-axisymmetric p-modes and gf-modes with n > 

^ . are damped and stable even in an inviscid disk. This eliminates a promising 

■ explanation for quasi-periodic oscillations in neutron-star and black-hole X-ray 
binaries. 

Subject headings: accretion, accretion disks — instabilities — waves — X-rays: 
binaries 



1. Introduction 

In recent years, the study of quasi-periodic oscillations (QPOs) in X-ray binaries has 
developed into a major field. A variety of QPOs have been observed in the variability power 
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spectra of neutron-star and black-hole X-ray binaries (van der Klis 2000; Remillard et al. 
2002, and references therein), and the observations have revealed a rich phenomenology. 

Of particular interest are the "kilohertz QPOs," which have frequencies of several hun- 
dreds of Hz to occasionally more than 10^ Hz. Because of their high frequencies, these QPOs 
must be produced by processes close to the accreting mass. However, since they have been 
seen in both neutron-star and black-hole X-ray binaries, it appears that the oscillations are 
not associated with the surface of the accreting object. Instead, it is generally believed that 
the kilohertz QPOs originate in the accretion flow surrounding the central mass. 

A detailed understanding of the oscillation modes of accretion disks could in principle 
allow observations of QPOs to be used to test strong gravity in the vicinity of compact 
objects (e.g., Stella & Vietri 1998; Stella, Vietri & Morsink 1999). Observations could also 
be "inverted" to measure relativistic parameters of the accreting mass, such as the spin of 
a black hole (Nowak et al. 1997; Wagoner, Silbergleit & Ortega- Rodriguez 2001). But such 
applications require a robust method of associating different observed QPO frequencies with 
specific disk modes and of calculating the frequencies of those modes from first principles. 

In what follows, we focus on gf-modes and p-modes in hydrodynamic thin disks. Follow- 
ing Silbergleit, Wagoner, & Ortega- Rodriguez (2001), Kato (2002), and references therein, 
we distinguish between the two kinds of modes primarily in terms of where most of the wave 
action is concentrated. If the bulk of the action is near corotation, we call it a gi-mode, and 
if the action is mostly away from corotation, we call it a p-mode. 

In an important paper, Okazaki, Kato, & Fukue (1987) showed that axisymmetric g- 
modes are trapped in the inner regions of a relativistic disk, where the epicychc frequency 
K reaches a maximum. This idea was exploited by Nowak & Wagoner (1991, 1992) and a 
number of other workers (Perez et al. 1997; Silbergleit et al. 2001; Abramowicz & Kluzniak 
2001) who worked out the physical properties of these and related modes (see Kato, Fukue 
& Mineshige 1998, Wagoner 1999, Kato 2001a for reviews). The modes are not dynamically 
unstable, however. Ortega-Rodriguez & Wagoner (2000) claim that viscosity destabilizes p- 
and g- modes, but we limit ourselves to inviscid disks. 

Recently, Kato (2001b, 2002) claimed to demonstrate that non- axisymmetric g-modes 
are trapped between two forbidden zones that lie on either side of the corotation radius. 
Kato further argued that these modes are highly unstable. This is very interesting because 
the observed kilohertz QPOs often have quite a large amplitude of flux variations, suggesting 
that the corresponding disk modes are probably dynamically unstable. This makes any disk 
modes that are dynamically unstable to be of great interest for the QPO problem. Although 
in a later paper, Kato (2003) has withdrawn the claim of an instability, his work nevertheless 
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suggests that it might be fruitful to explore dynamical instabilities in non-axisymmetric disk 
modes. This is the motivation behind the present paper. 

The work presented here is particularly influenced by the papers of Goldreich & Narayan 
(1985, henceforth GN) and Narayan, Goldreich, & Goodman (1987, henceforth NGG). These 
authors studied perturbations in a particular simplified fluid system called the shearing sheet. 

They showed that non-axisymmetric modes with no vertical nodes (n = in the notation 
used in this paper) can be dynamically unstable. They further demonstrated that the driving 
force for the instability is a wave amplifier inside the system. These papers were in turn 
influenced by work on spiral wave amplification in stellar disks, especially Mark (1976). 

We use the physical understanding obtained from the above work to guide our present 
investigation. In §2 we consider a thin accretion disk (which is more general than the 
shearing sheet), and derive a wave equation for linear non-axisymmetric perturbations and 
an associated conserved current. In §3 we explore the nature of the singularities in the wave 
equation associated with the corotation resonance and the Lindblad resonances. In §4 we 
consider the WKB limit of the wave equation and identify the basic properties of ingoing and 
outgoing waves in various regions of the disk. In §5 we consider the interaction of waves with 
various barriers associated with the Lindblad and/or corotation resonance. We show that 
waves with n = are amplified when they reflect off the corotation barrier. This confirms 
the result obtained by GN and NGG. However, when we repeat the analysis for waves with 
n > (waves with one or more nodes in the vertical direction), we find that there is no 
amplifier anywhere in the system, either at the Lindblad resonances (§5) or at corotation 
(§6). Indeed, we show that corotation acts as a severe absorber of waves. Based on these 
results, we conclude in §7 that there are no dynamically unstable p- or g^-modes with n > in 
a thin disk. Appendix A presents some numerical results which help to extend the analysis 
to the non-WKB regime. 

2. Wave Equation for Non-Axisymmetric Perturbations to a Thin Disk 

We consider an axisymmetric thin disk. Because of the symmetry of the unperturbed 
system, we assume without loss of generality that the linear perturbations are proportional 
to exp[i{—u!t + m0)] in cylindrical coordinates (r, (f), z). Here, m is an integer and uj is the 
mode frequency, which may be either real or complex. We write the linear perturbations to 
the mass density and velocity as 



p{t,r,(f),z) 




Po{r, z) + pi(r, z) exp[i(-a;t + m0)] , 




(2) 



v^{t, r, (f), z) 
Vz{t, r, 0, z) 



^{r)r + v{r, z) exp[i(— + m0)] , 
w{r, z) ex^[i{—ujt + m^)] , 



(3) 
(4) 



where po is the unperturbed mass density, Vt is the unperturbed angular velocity, and pi, 
M, w, and w are first-order perturbations to the mass density and the velocity field. Note 
that dfl/dz = since we take the disk to be isentropic. Following Kato, we use Newtonian 
equations. Although general relativity may introduce additional instabilities, many of its 
effects can be mimicked by appropriate radial variations of the structural disk frequencies 
Q, K, and Q_L, the latter two of which are defined below. 

For simphcity, we assume that the disk has an isothermal equation of state. 



where p is the gas pressure and Cg is the sound speed which we take to be a constant. The 
equilibrium structure of a sufficiently thin isothermal disk then takes the form 



where h — h{r) — Cs/fl±{r) is the half-thickness of the disk, and Q± is the vertical disk 
frequency, which may be different from the orbital angular velocity due to strong radial 
pressure gradients in a Newtonian disk (unlikely for a thin disk) , or to relativistic gravity. 

The isothermal equation of state allows an almost exact separation of variables, as 
shown below. This simplification is well worth the loss of generality, especially since the 
main focus of our paper is on the corotation resonance. In the vicinity of the corotation 
resonance, in fact, the frequency a (defined below) is small compared to the reciprocal 
of the sound crossing time across the thickness of the disk; consequently, motions near 
corotation are noncompressive, so that the choice among different isentropic equations of 
state is not important. Far from corotation (|r — rd > r/m), however, non-isothermal 
equilibria concentrate waves toward the disk surface where they may be more easily dissipated 
(Bate et al. 2002, and references therein). 

In terms of the enthalpy perturbation 



(5) 




(6) 




(7) 



the first-order perturbation to the continuity equation is 
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where 



a = uj — mfl 



(9) 



is the frequency of the perturbation in the local corotating frame of the disk. The first-order 
perturbations to the radial, azimuthal, and vertical momentum equations lead to 



respectively, where 



dQ 

iau + 2uv — —— , 
or 

■^u-iav = Q , 

2\L r 



law — 



dQ 
dz 



(10) 

(11) 
(12) 



dr 



(13) 



is the epicyclic frequency of the flow in the r0-plane. 

Solving for u and v from equations (10) and (11), and w from equation (12), and sub- 
stituting the results into equation (8), we obtain a second-order partial differential equation 
for Q 



a"^ d f por dQ 



Por dr 
where 



D dr 



cr^ ^ ^ma\^^ 1 2ma d 



\ r J D ' Por dr \ D 



+ 



D = i<^-a'' 



(15) 



Equation (14) is most easily analysed by separation of variables in r, z, as discussed by 
Kato (2001a). Since dh/dr ~ h/r is small for a thin disk, we neglect the dependence of h on 
r and introduce a coordinate rj = z/ [V2hj to replace z. Then, writing Q{r, rj) = Qri^)Qriiv) 
and substituting in equation (14), we can show that Qr{r) and Qniv) satisfy 



d_ f Poor dQ, 
dr \ D dr 



Poor 



a^ 



m'^a'^ 



2ma d 



f Poo^ \ 
r'^D Poor dr \ D J 



= , 



(16) 



where n is a constant. 



(17) 
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Equation (17) has the form of the Hermite equation. In order for the energy density 
of perturbations to be bound as r) ^ ±00, n must be a non-negative integer (Okazaki et 
al. 1987): n = 0,1,2,.... Then the solutions to equation (17) are given by the Hermite 
polynomials: = Hn{ri). The integer n determines the number of nodes in the vertical 
direction. An even (odd) n corresponds to an even (odd) mode of oscillation. The mode 
with n = has Qn = constant, and has no motions in the vertical direction according to 
equation (12). 

Equation (16), with n = 0, 1,2, is the wave equation for linear perturbations of an 
isothermal thin disk. Kato (2001a) obtained a similar equation for more general equations 
of state by neglecting slowly radially varying terms. 



The quantity 



2W \ dr dr J poor 



where the asterisks denote complex conjugates, represents a current of wave action. It can 
be checked that this current is conserved, that is, dJ/dr — 0, wherever equation (16) is 
nonsingular and uj is real. The current plays a prominent role in all of the analysis presented 
below. The conservation law for wave action can be extended to complex frequencies uj — 
ujR + iuji, 

dJ 
dr 

where pa is a real quantity proportional to |QrP representing the density of wave action 
(NGG). The role of the singularities is more important in the present context. We shall 
show that J is conserved at all singularities except corotation, and even there it is conserved 
when n — 0. 



3. Corotation and Lindblad Singulatrities in the Wave Equation 

When a is real (i.e., cu is real), the wave equation (16) contains two types of singularities: 
one is the corotation singularity given by the condition a — 0, the other is the Lindblad 
singularities given by the condition D — 0. The former occurs at the corotation radius Tc 
and the latter at the Lindblad radii tl. 
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3.1. Corotation singuleirity: n — 

We begin by first discussing the case n = 0, which corresponds to the problem analysed 
by GN and NGG. They considered the special case of the shearing sheet, whereas we consider 
here a more general disk; however, near corotation, the two problems are very similar. For 
r close to Tc, the wave equation becomes 



(PQr 



TcXr - Tc 



0, 



— dlnO/dlnr dlnr \ D J 



(19) 



where we have used a = m {—dfl/dr)^^^^ {r — rc) near r = Tc ■ Assuming that Qr oc (r — Tc)^^ 
near r — rc, the leading term in equation (19) gives rise to 



/3i(/3i-l)(r-rc)'^i-2 = 



(20) 



which requires /3i = 0, or 1. Thus, near the corotation radius the solutions to the wave 
equation are 



oc 1 , 



(21) 



Equation (21) implies that the solutions are analj^ic near r = Tc- Correspondingly, as 
can be easily checked, the current defined by equation (18) is conserved across Tc, i.e., 

^ = 1 , J^^ = J (r - re = 0±) . (22) 

Therefore, when n — 0, there is no real singularity at corotation, only an apparent singularity. 
Indeed, for the model of a shearing sheet, NGG have shown that the corotation singularity 
does not appear at all if the perturbation equation is written in terms of the azimuthal 
velocity rather than Qr- 



3.2. Corotation singuleirity: n > 



When n > 0, the wave equation near the corotation radius becomes 

Q I K 



d'^Qr _^ nb'^Qr 



0, 



dr^ (r — Tc)^ 
The solutions to equation (23) are 

Qr {r - rc)-^^^'\ 



b = 



mcs{—dn/dr) 



(23) 



r=rc 



(24) 
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Due to the presence of the factor (r — r^)^/"^, the solutions are not analytic at r = Tc- Thus, 
we have a real singularity. Note that, for a thin disk, we usually have h ~ r/mh ^ 1 unless 
m is > r//i , so generally we expect g ^ 1. 

Substituting the two solutions in (24) into equation (18), we can calculate the cor- 
responding current density Jc+ defined in the previous subsection. Making the analytic 
contimiation r — — > (r — rc)e'^' to the solutions in equation (24) and substituting the 
results into equation (18), we can calculate the corresponding current Jc_. We find that the 
current is not conserved across corotation: 



Jc- 



c+ 



= e^^'^^ , (25) 



where the upper/lower sign corresponds to the upper/lower sign in equation (24), respec- 
tively. 

This analysis shows that the behavior near corotation is very different for n = and 
n > 0. This is one of two major differences between the two cases, the other being the 
geometry of permitted and forbidden zones for waves, discussed in §4 and displayed in 
Figures 1 and 2. 



3.3. Lindblad singuleirity 

The Lindblad singularity occurs at the radii tl where D — Q. Near this singularity, we 
have a ^ ±k and D D'i^{r — tl), assuming that D'l^ = {dD / dr)r=ri^ 7^ . Then, assuming 
that Qr oa {r — ri^Y"^ near r = tl, the leading term in equation (16) gives rise to 

/?2(/92-2)(r-rL)^^-3 = 0, (26) 

which requires /32 = 0, or 2. Thus, when Dl 7^ , the two solutions to the wave equation are 
locally 

a 1 , (r - r^f . (27) 

Both solutions are clearly analytic. 

More generally, if D is an analytic function of r near the Lindblad radius (as it is in our 
problem), the solutions to the wave equation must also be analytic there. It can be checked 
that the current defined by equation (18) is then conserved across the Lindblad singularity. 
Therefore, the Lindblad resonance is not a true singularity but only an apparent singularity 
(NGG; Kato 2002). This is self-evident when the second-order equation (19) is replaced by 
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the equivalent system of first-order equations (A12). Indeed, if one considers wave equations 
for other perturbed quantities (e.g., radial velocity), the Lindblad singularity gives way to 
other singular terms, which occur at different radii and are again not real singularities but 
only apparent. The corotation singularity, on the other hand, persists for any fluid variable 
one selects and is a genuine singularity when n > 0. 

In summary, when n > the corotation resonance is an intrinsic singularity in the 
wave equation, where the conservation of the current breaks down. When n = 0, the 
corotation resonance is only an apparent singularity; the solutions are well-behaved there 
and the current is conserved across the resonance. The Lindblad resonance is never a real 
singularity; the solutions are well-behaved and the current is conserved across the Lindblad 
resonance for all values of n. 



4. WKB Solutions to the Wave Equation in Permitted Regions 



In the WKB regime, where the wavelength is much smaller than the length scale over 
which the potential in the wave equation varies, the wave equation (16) can be approximated 

by 



fOr _ dlnWdQr 
dr"^ 



dr dr 



+ k^,{r)Qr^O, 



where 



(28) 



(29) 



In the last step, we have assumed that mh/r ^ 1. Equation (28) is Kato's (2001b, 2002) 
wave equation, except that we have retained the subdominant term involving the Wronskian 
W (eq. [18]). This term is crucial because it influences the sign of the wave action and hence 
the amplification of waves. 



When the condition 



1 dk. 




1 


kl dr 




kj-r 



< 1 



(30) 



is satisfied, the wave equation (28) can be solved with the WKB approximation (see, e.g., 
Merzbacher 1998): 



Qr, 



w 



WKB 



exp ( ±i / krdr 



(31) 
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The solution with the "+" sign in the exponential corresponds to wave- vector -l-A;^, and the 
solution with the "— " sign corresponds to wave-vector —kr- 

Within the WKB approximation, permitted regions for waves are defined by the re- 
quirement k^{r) > 0, and forbidden regions by k^{r) < 0. Consider first the simpler case of 
n — 0. Here, the region of the disk in which cr^ < is forbidden. This region extends from 
the inner Lindblad radius {a = —k) to the outer Lindblad radius {a = k), straddling the 
corotation radius ((7 = 0). It is shown as forbidden region A in Figure 1. On either side of 
A, there are two permitted regions: region I with a < —k, and region II with a > k. 

When n > 0, things are more complicated. Now, there are four permitted regions (see 
Fig. 2): 

- Region I: the region with a < — max {^/nfl±, n) ; 

- Region II: the region with a > max {y/nft±, k) ; 

- Region III: the region with — min {^/nQ±, k) < a < 0; 

- Region IV: the region with < cr < min {^/nfl±, k) . 

Regions I and III have a < and are on the left-hand side of corotation. These two zones 
are separated from each other by a forbidden zone (barrier Al, see Fig. 2). Regions II and 
IV have a > and are on the right-hand side of corotation. They are again separated from 
each other by a forbidden zone (barrier A2). Regions III and IV are separated from each 
other by the corotation resonance, which is not a true barrier. In contrast to the n — case, 
here the region around corotation is permitted. 

For a WKB solution with a wave- vector kr (i.e., the solution in eq. [31] with the 
sign in the exponential), the corresponding current is 

kr 2 

■AvKB ~ |Qr,WKB| , (32) 

where we have used the fact that in the permitted regions kr is real so that k* — kr- From 
the definition of W (eq. [18]), the sign of W is determined by the sign oi D — — a^. Hence, 
the current and the wave- vector have the same signs in regions where > (regions I and 
II, see Table 1), but opposite signs where cr^ < k,^ (regions HI and IV). 

We now define an "outgoing" wave as one that moves away from corotation, i.e., towards 
larger \a\, and an "ingoing" wave as one that moves towards corotation (see Figs. 1 and 2). 
The direction of motion is determined by the sign of the group velocity Vgr rather than the 
wave- vector kr- 
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By definition, Vgr = {dkr/duj) — (dkr/da) , so that 



din fe 



-1 



^rVgr = (T ( ) . (33) 



Prom equation (29) we have 



(34) 



9 In (72 ((72 - ) ((72 _ ^2) ' 

and therefore 

= (7^ - nniK^ ■ ^^^^ 

It is then easily verified that kj.Vgr > in regions II and III, and /ovfgr < in regions I 
and IV. This allows us to identify whether a particular wave in a given permitted region is 
ingoing or outgoing. 

Knowing the signs of krJ and krVgr in each permitted region, we are able to determine 
the sign of VgrJ in that region (Table 1). This is an important quantity, as we show in the 
next section. 

Finally, we briefly comment on the locations of the four permitted regions in a real disk. 
From their definitions, if all four regions I, II, III, and IV exist, they should appear in the 
following order with increasing disk radius (as in Fig. 2): region I, region III, region IV, and 
region II, since a is an increasing function of radius. Thus we expect region I to be near the 
inner boundary of the disk (r ~ Tin), region II to be at large radii (r ^ Tc), and regions III 
and IV to lie in between, on either side of corotation. 



5. Wave Ampliflcation/De-amplification at the Bcirriers 

As explained in §1, a dynamical instability is possible if there is a wave amplifier in 
the system. In the case of n = 0, such an amplifier is present at the corotation barrier 
(GN, NGG), as we now demonstrate using the results derived in §4. Figure 1 shows that, 
when n — 0, there are only two permitted regions in the disk: region I near the inner 
boundary of the disk, and region II at large radii. These two regions are separated from 
each other by barrier A, which includes the corotation and Lindblad resonances. However, 
as shown in §§3.1 and 3.3, neither the corotation resonance nor the Lindblad resonances are 
real singularities, and the current defined by equation (18) is conserved throughout the disk. 

Consider now a wave incident on the potential barrier from region I. The interaction of 
this ingoing wave with the potential barrier will produce a reflected outgoing wave in region 
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I and a transmitted outgoing wave in region II. Because of the conservation of the current, 
we have the following relation among the currents of the three waves: 



'-'ir, ' "-'f 



out 



x 



(11) 



out 



(36) 



From Table 1, in region II the group velocity and the current density of the wave have 
the same sign. Since an outgoing wave in Region II has a positive group velocity (the wave 
goes away from corotation towards larger r) , the corresponding current J^^^ must be positive. 
Thus, the right-hand side of equation (36) is positive. In region I, the group velocity and the 
current of the wave have opposite signs. The ingoing wave here has a positive group velocity 
(it is moving from small values of r towards Tc), while the outgoing wave has a negative 



group velocity. Thus, J^^^ < and J^^^ > 0. Then, equation (36) is equivalent to 



J, 



(I) 



out 



(I) 



J, 



(11) 



out 



(37) 



We define the gain G to be the absolute value of the reflection coefficient. Equation (37) 
implies that the gain 



G 



Jr 



reflected 



J incident 



J, 



(I) 



(I) 



>1, 



71 = 0. 



(38) 



In other words, for n = 0, an incident wave is reflected by the potential barrier with a larger 
amplitude, and the potential barrier behaves as an amplifier. If, in addition, the inner edge 
of the disk behaves like a near-perfect refiector, then we have the makings of an instability. 
(This would be a p-mode disk instability since the wave action is concentrated away from 
corotation.) 

The corotation amplifier works equally well if a wave is incident on the potential barrier 
from region II, as is easily verified. In this case, to have an instability, the outgoing wave 
in region II must be somehow refiected. The refiection is unhkely to be from the distant 
outer edge of the disk. Perhaps an inhomogeneity in the disk might provide the necessary 
reflection, but the topic is beyond the scope of this paper. For the purposes of this paper, 
the key point is that, when n = 0, the system has a wave amplifier and, therefore, can 
potentially have unstable modes. 

The situation changes dramatically when n > 0. Now, there are four permitted regions 
in the disk, 1, II, III, IV, and two barriers, Al, A2. Let us consider the region to the right 
of corotation (i.e., o" > 0). Assume that an outgoing wave is incident on the barrier A2 from 
region IV. The interaction of this outgoing wave with the barrier produces a reflected ingoing 
wave in region IV and a transmitted outgoing wave in region II. We have shown earlier that 
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there is no real singularity at the barrier, and so the current density must be conserved from 
region IV to region II. We then have 



j(IV) J 



(IV) 



J, 



(11) 



out 



(39) 



Prom Table 1 we see that in region IV the group velocity and the current have the same 
sign. Since in region IV an outgoing wave has a positive group velocity and an ingoing wave 
has a negative group velocity, we have J^^^ > and Jj^^^ < 0. Similarly, J^"-* > 0. Thus, 
equation (39) can be rewritten as 



J 



(IV) 



7(IV) 

■^out 



J, 



(11) 



out 



(40) 



Equation (40) implies that the gain in this case is always less than unity: 



G 



J reflecti 



cd 



J incident 



J; 



(IV) 



J. 



(IV) 



out 



< 1, 



n > 0. 



(41) 



In other words, a wave that is incident on barrier A2 from region IV is always reflected with 
a smaller amplitude. It can be easily checked that this statement is true also for a wave 
incident on barrier A2 from region II, and also for waves incident on barrier Al from either 
region I or region III. 

We thus conclude that, for n > 0, neither of the two barrier Al and A2 behaves hke 
an amplifier; both barriers deamplify waves. This eliminates a promising mechanism for 
producing a dynamical non-axisymmetric instability. The only thing left to be checked is 
whether corotation itself can amplify waves. This is the topic of the next section. 



6. Absorption at the Corotation Resonance 

To understand the role played by the corotation singularity when n > 0, we need to 
study the behavior of the solutions of the wave equation (16) near a — Q (i.e., r = r^)- 
In particular, we would like to know for a wave incident on the resonance, how much is 
transmitted to the other side and how much is reflected. 

As |(t| — > (i.e., r — > Tc), and taking n > 0, the wave equation becomes equation (23), 
whose two linearly independent solutions are given by equation (24) for q ^ 0. The wave- 
vectors corresponding to the two solutions are 



(42) 
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It can be checked that when ^Jnh » 1 (i.e., ? ^ 1) the solutions in equation (24) become 
the WKB solutions given by equation (31). 

According to Table 1, in region IV (r > the wave- vector and the group velocity have 
opposite signs. Thus, in region IV the ingoing wave, which has a negative group velocity 
i.e. a positive wave- vector, is given by (r — rc)2+^**. Now let us analj^ically continue this 
sohition into region III. As discussed in NGG, the continuation must be done such that the 
integration path in the complex plane goes above the singularity at a = if the solution is 
to correspond to an initial-value problem. That is, we must take r — ^ \t — rd e'^*. We 
then find that the solution transforms as follows: 

Region IV ^ Region III : (r - rc)^+^''' ^ ie'^'^^ |r - rd^^^^'^ . (43) 

In region III (r < the wave-vector and the group velocity have the same sign. So, in 
region III the outgoing wave, which has a negative group velocity i.e. a negative wave-vector, 
is given by {r — r^^^'^'^'^ oc |r — rd^^^**. Thus, from equation (43), the ingoing wave in region 
IV becomes a purely outgoing wave with a reduced amplitude in region III, and there is no 
reflected wave. 

Similarly, for an ingoing wave in region III we find 

Region III ^ Region IV : {r^ - r)^-^''^ ^ -ie'^^^ |r - rcp"^'* , (44) 

where we have taken Tc — r — > |r — rd e""^'. 

Computing currents, we may calculate the transmission and reflection coefficients T^. 
and Rc of the corotation resonance: 

Te = e-^"^ i?e = . (45) 

Thus, there is no wave reflection at the corotation singularity, and there is a severe absorption 
of wave action in the transmitted wave. Indeed, the absorption is exponentially strong when 
q is large (as for a thin disk). 

From equation (45), it might appear that the absorption would disappear for nl? < 1/4 
since then q becomes zero or imaginary and so |Td = 1. However, this is not true in general, 
and absorption disappears only when nl? = (i.e., q = i-i/4:). For example, consider a 
solution in region IV containing both ingoing and outgoing waves 

Q(iv) ^ A,{r - r,)-^+^"i + A2{r - re)^-^^^ , (46) 

where Ai and A2 are complex numbers. The corresponding solution in region III obtained 
by analytic continuation is 

Q("^) = iAie-''^^ |r - r,\ ^+''^ + ^ae'"^ |r - r,\ ^""'^ . (47) 
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Then we can calculate the net currents on the two sides of corotation. We find 

.(IV) .(IV) , 
•^in ~r ^out i 

_(- j^™) COS Anq + cot V' sin Anq ' 

where q = iq and 2-?/' = arg(AiA2/A^A2). The ratio of the currents is equal to 1 for any 
if and only if q = ±i/4:, i.e. nb^ = (in which case the corotation singularity disappears 
or becomes an apparent singularity and we return to the case studied by NGG). Thus we 
conclude that, in general, non-axisymmetric gf-modes in disks are absorbed at corotation. In 
other words, the corotation singularity de-amplifies waves and therefore cannot induce an 
instability. 

Absorption at corotation is not unique to disk gi-modes. It occurs for other sorts of 
waves in shear flows and has been well studied in the fluid-dynamics literature, e.g. Drazin 
& Reid (1981). In particular, a nonrotating, incompressible, stratified shear fiow is absolutely 
linearly stable when the Richardson number 

_ —gd\np/dx 
" {dVy/dxy 

is greater than 1/4 (Howard 1961). Here g is the acceleration of gravity, Vy the unperturbed 
velocity, and x the vertical direction. Booker & Bretherton (1967) studied gf-modes in such 
flows as an initial value problem and showed that all but a fraction exp (^—ny^R — 1/4^ of the 
wave current is absorbed near the altitude Xc where V{xc) = ur/ky, called the "critical layer". 
They also showed that in the WKB approximation, wave packets propagate towards the 
critical layer but fail to reach it in flnite time. In the present problem, the role of the critical 
layer is played by corotation, and that of the Richardson number by nb"^ — uk^/ (mhdQ./dr)'^ . 

Combined with the results in §5, we see that, when n > 0, there are no amplifiers in the 
system, either at the two barriers or at corotation, and so there are no non-axisymmetric 
growing modes trapped in the disk. This result has been obtained here via a WKB analy- 
sis. We confirm the result in Appendix A with numerical calculations for the specific case 
of the shearing sheet model. That analysis is more general and goes beyond the WKB 
approximation. 



(48) 



7. Summary and Discussion 

We have derived the wave equation for non-axisymmetric perturbations in an isothermal 
hydrodynamic thin disk, and have identified a conserved current (§2). The wave equation 
contains two types of singularities: the corotation singularity and the Lindblad singularity. 
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The Lindblad singularity is never a real singularity, regardless of the vertical wave-number 
n; the solutions are always analytic and well-behaved in the vicinity of this singularity, and 
the current is conserved across it (§3). For n = (horizontal motions independent of height), 
corotation is also not a real singularity, and the current is conserved. But for waves having 
vertical nodes {n > 1), corotation is a true singularity where the conservation of current 
breaks down. 

Using a WKB approach (§4), we have shown that for n = there are two permitted 
regions in the disk (Fig. 1), one near the inner boundary (region 1), and the other at large radii 
(region II). In contrast, for the first and higher overtones (n > 0), there are four permitted 
regions in the disk (Fig. 2): region I near the inner boundary, region II at large radii, and 
regions III and IV in between, on either side of corotation. We have analyzed the WKB 
solutions in each permitted region, mapped the ingoing and outgoing waves, and identified 
the signs of the current density and group velocity of the various waves (summarized in 
Table 1). 

For n — 0, the current is conserved throughout the entire disk, from the inner boundary 
to infinity. This, combined with the results given in Table 1, allows us to demonstrate that 
the disk behaves like a wave amplifier (§5). The argument is simple. We consider an ingoing 

wave in region I, which produces a reflected outgoing wave in region I and a transmitted 
outgoing wave in region II. Since the outgoing wave in region II has a positive current, by 
the conservation of current, this means that the sum of the ingoing and outgoing waves in 
region I should also have a positive current. However, in region I, the outgoing wave has a 
positive current and the ingoing wave has a negative current. Thus, the outgoing wave in 
region I must have a larger amplitude than the ingoing wave. In other words, the interaction 
of the ingoing wave with the corotation region has caused wave amplification in the reflected 
outgoing wave. 

This result for n = waves was shown for the case of the shearing sheet model by GN 
and NGG. Indeed, these authors showed that the amplifier causes an instability if a reflecting 
boundary condition exists in one of the permitted regions. The reflective boundary causes 
the amplifled outgoing wave to return to corotation and be amplified repeatedly. While this 
instability is a candidate to explain QPOs in accretion systems, it unfortunately requires 
nearly perfect reflection at the boundary, presumably the inner edge of the disk, which is 
somewhat problematic. As the analytical results in the above-quoted papers indicate, and 
also confirmed in the numerical calculations of Appendix A of the present paper (see the first 
two rows of Table 2), the amplification at the corotation barrier is usually extremely weak 
for a thin disk. Thus, any energy loss either during transit of the wave to the boundary and 
back, or during the reflection at the boundary, would kill the instability. The requirement of 
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perfect reflection at the boundary is particularly troublesome, since Blaes (1987) has shown 
that any radial inflow of the gas at the boundary (which is unavoidable in an accretion flow) 
would severely reduce the reflectivity of the boundary. 

What is needed is a sufficiently strong amphfler, such that even less than perfect reflec- 
tion at the boundary is sufficient to give overall amplification and instabihty. Kato's (2001b, 
2002) claim that non-axisymmctric (^-modes with n > are strongly amplified appeared 
to be just the answer. Unfortunately, Kato (2003) himself discovered that (jf-modes do not 
grow, and our independent analysis described in this paper confirms his result. We have 
shown in §4 that for first and higher overtone modes (n > 0), there are two barriers, one 
near each Lindblad resonance, instead of the single barrier for the n = case (compare Figs. 
1 and 2). Both barriers unfortunately behave as wave de-amplifiers rather than amplifiers 
(§5); for a wave that is incident on either barrier from either side, both the refiected wave 
and the transmitted wave are weaker than the original wave (the total current is however 
conserved). Furthermore, for these n > waves, the corotation singularity behaves as a 
strong wave absorber (§6) so that current is absorbed and is lost whenever a wave is incident 
on corotation. Thus, we conclude that n > modes are unable to grow in a thin disk. 

If we define gi-modes as those that have the bulk of their wave action near corotation 
(Silbergleit et al. 2001; Kato 2002), then these modes must of necessity have n > 0, since 
only such modes have a permitted region near corotation (Fig. 2). Our analysis thus shows 
that non-axisymmetric (?-modes cannot be dynamically unstable, p-modes, which have their 
wave action far away from corotation, are possible both for n = and n > 0. Unstable 
n > p-modes are again ruled out by our analysis. This leaves only the n = modes, which 
we discussed earher, as candidates for QPOs. We have already argued that these modes 
are unlikely to grow in thin disks. Perhaps with a sufficiently thick disk, the wave amplifier 
could become strong (e.g., see the third row in Table 2) and might give an instability. But 
this is not very promising since it is likely that a thick disk will also have more rapid gas 
infiow at the inner boundary and therefore weaker refiection there. 

Perhaps by adding more physics to the disk model one may yet find a global hydrody- 
namic instability suited to explain QPO. As already noted, axisymmetric viscous instabili- 
ties have been proposed by Ortega-Rodriguez & Wagoner (2000). Vertical stratification (i.e., 
perpendicular to the mean fiow) may allow non-axisymmetric instability in some laboratory 
Couette flows (Yavneh, McWilliams, & Molemaker 2001). In our opinion, neither of these 
is a very plausible instability mechanism for disks, but space does not permit an adequate 
discussion here. 

The answer may lie with magnetohydrodynamic effects. To date, most detailed studies 
of magnetorotational instability have emphasized the production of turbulence rather than 



-18- 



high-Q global oscillations, but we suspect that this is the most promising direction for future 
work. 
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A. The Shearing Sheet Model 

In this Appendix we approximate a differentially rotating thin disk by means of the 
shearing sheet, which consists of a uniform shear flow with a Coriolis force (NGG). We use 
a Cartesian coordinate system {x, y, z), where x and y are related to the polar coordinates r 
and (j) of the original disk by 

X = r - Tc , y = rc(0 - ftct) ■ (Al) 

Here, Tc is a reference radius and flc = ^(^c) is the disk angular velocity at Tc. The velocity of 
the unperturbed flow is v = 2Axj where j is a unit vector in the y-direction. The frequency 
A and the related frequency B are the Oort constants of the disk at Vc- 

(A2) 

r=rc 

The epicyclic frequency of the flow in the xy-plane is 

K = ^/ABn^ = y/AB{B - A) . (A3) 

We treat the frequencies A, B, k as constants. 

When Q oc where g is a constant, we have 2A — —qi^c: 2B — {2 — q)flc: a-nd 
K — •\/2(2 — q) fic- It is useful to remark that q — 2 corresponds to a disk with constant 
angular momentum, g = 3/2 to a thin Keplerian disk, and g = 1 to a disk with constant 
circular velocity. For < g < 2, we have A < and B > 0. 

We assume that the flow is isothermal with an equation of state p — pc^, where p is the 
gas pressure, p is the mass density, and Cg = constant is the sound speed. We take the flow to 



A = 



rdn 

2~d^ 



2r dr 
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be homogeneous and to extend to ±00 along x and y. In the vertical direction we make the 
usual harmonic approximation for the potential, with vertical frequency Qj^. Then, vertical 
hydrostatic equilibrium implies p (x p (x exp(— where the vertical scale height h is 
equal to Cs/Q±. In a thin Keplerian disk, fl± — flc — but we allow these three frequencies 
to be different. 

Without loss of generality, we assume linear perturbations proportional to exp[i(— + 
kyy)], where ky = m/rc is the azimuthal wave- vector, and the frequency cu may in general 
be complex. Then, the first-order perturbation equations are 

QQ 

iau + 2D,cV — , (A5) 
2Bu — iav — —ikyQ , (A6) 

= ' (A7) 

where u, v, and w are the perturbations to the velocities in the x, y, and z directions 
respectively, po is the unperturbed mass density, Q = Sp/po is the perturbed enthalpy, and 



UJ 



mn{r) = -2kyAx (A8) 



is the pattern frequency of the mode, where we have chosen Tc to be the corotation radius 
where the fluid moves with the same speed as the pattern speed of the mode 

mO(rc) = UJ . (A9) 

Equations (A4-A7) correspond to equations (8) and (10-12) for the case of a thin disk, 
respectively. 

Prom equations (A4-A7) a second-order partial differential equation in x and z may 
be derived for any of the perturbed flow quantities. Following Kato (2002), we consider 

the equation for Q. The isothermal equation of state allows a separation of variables, Q = 
Qx{x)Hn{z/ \/2h), where Hn{i]) is a Hermite polynomial of degree n > (sec Okazaki et al. 
1987 and Kato 2002 for details). The function Qx{x) then satisfies the following differential 
equation: 

where D = — a^. Equation (AlO) is similar to equation (16) which we derived for the 
thin disk. In the same way, we can define a conserved current by 
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As in a thin disk, there are two types of singularities in equation (AlO): the Lindblad 
singularity at D = 0, and the corotation singularity at a = 0. The Lindblad singularity 
is always an apparent singularity. The corotation singularity is a true singularity when 
n > but apparent when n = 0. In the WKB limit, there are four permitted regions 
in a shearing sheet flow when n > 0: region 1 with a < —ma.x{^/n^l±,K), region II with 
a > max{^/nfl±, k), region III with — min(Y^^^±, n) < a < 0, and region IV with < cr < 
mm{^/nQ±, k.). Region I (II) and region III (IV) are separated from each other by a potential 
barrier, region III and region IV are separated from each other by the corotation singularity. 
When n = 0, there are only two permitted regions: I and II. 

The results in the main text can be applied to the shearing sheet model. In particular, 
when n — 0, the potential barrier between regions I and II behaves as an amplifier, a wave 
incident on it is reflected with larger amphtude (NGG). When n > 0, the two potential 
barriers (one between regions I and III, the other between regions II and IV) behave as a 
de-amplifler, a wave incident on it is reflected with smaller amplitude. In addition, when 
n > 0, the corotation singularity absorbs energy from a wave. 

Here we show some numerical results for the shearing sheet model to confirm these 
results. Since these numerical solutions do not require that the WKB approximation be 
valid in regions III and IV, they complement the results obtained in the main text. 

For numerical work, it is convenient to replace equation (AlO) with an equivalent set of 
two first-order differential equations. Assuming as before that all variables depend upon z as 
Hermite polynomials Hn{z/\/2h) or its derivative, we eliminate w and v from equations (A4)- 
(A7) to obtain, in matrix form. 



d 




1 


-2kyB 






Ux 


dx 


Qx 


a 


_ i{a^ - «2) 


2 Jvy^^ Q 




Qx 



(A12) 



In this first-order form of the linearized equations, it is clear that the only possible singularity 
is 0" = 0. This explains why the Lindblad resonances a = ±k are not true singularities of 
the equivalent second-order equations (AID) and (16). 

We look for solutions to equations (A12) that contain only outgoing waves in region II. 
To do so, we start from an outgoing wave at xq > in the WKB region II, which has a 
positive wave- vector (thus a positive group velocity), then integrate the equations along the 
real axis of +x toward the corotation. Near but before crossing the corotation, we deform 
the integration path into the complex plane of x to pass the corotation singularity from 
above, then come back to the real axis of —x. Then, we integrate the equations along the 
real axis of —x into the WKB region I, until x — —xq is reached. At x — —xq, we decompose 
the solution into an ingoing component (having a negative wave- vector, i.e., a positive group 
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velocity) and an outgoing component (having a positive wave- vector, i.e., a negative group 
velocity). Then, we calculate the gain G (i.e., the reflection coefficient), and the transmission 
coefficient T by 



G 
T 



Ux,in{x = -Xq) 
Ux,out{x = Xo) 



Qx,out{x — —Xq) 



Ux,m{x = -Xq) 



Qx,in(^X — Xq) 
Qx,out{x — Xq) 



Qx,in(^X — Xq) 



(A13) 
(A14) 



where "out" denotes "outgoing", "in" denotes "ingoing". 



Numerical results corresponding to a few choices of parameters are shown in Table 2. 
They are classified into two classes: one with n = 0; the other with n = 1. Each class 
contains both a Keplerian disk (O^ = n = Q^) and non-Keplerian disks {Q± and k different 
from Qc)- The results show that, for the case of n = 0, where the corotation is not an 
intrinsic singularity, the incident wave is amplified: G > 1; indeed G — 1 + T since the 
current is conserved. One can check that the numerical results agree with the analytical 
results of NGG, where 



G = 1 + T = 1 + exp(-27rC) , 



C 



C^ky + 

4:Cs\Aky\ 



(A15) 



However, for the case of n = 1, where the corotation is an intrinsic singularity, there is strong 
absorption at the corotation as indicated by the fact that T -C 1 — G, and the incident wave 
is always de- amplified since G < 1. 

The initial and terminal points of our numerical integrations are always chosen well 
within the WKB domain of regions II and 1, respectively (i.e., \kxx\ ^ 1, where k^ is 
the wave- vector in the x-direction), so that the ingoing and outgoing waves are clearly 
distinguishable. It can be shown that for n > 1, the criterion for the validity of the WKB 
approximation in regions III and IV near corotation reduces to 6 ^ 1, where 



b = 



K 



-2Akyh 



The solutions listed in Table 2 do not satisfy 6^1. Indeed, for the last solution, where 
rt± = Omc, K = 0.2^]c, and kyh = 0.3, we have b = 0.336 and q = (1/2)^/^62 - 1/4 = 
0.185i. Even for this case with an imaginary q, we see that T < 1 — G, which must be 
attributed to absorption at corotation. These numerical results support and complement 
our analytical results based on the WKB approximation. 
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Table 1. Relations among the directions of the current, the wave- vector, and the group 

velocity in the four permitted regions 



Region D = — a"^ k^J^ k^v^r^ VgrJ^ 

I - + - - 

II - + + + 

III + - + - 

IV + - - + 



'^The sign of the product krJ determines the relative directions 
between the current J and the wave- vector kr- 

^The sign of the product krVgr determines the relative direc- 
tions between the wave- vector kr and the group velocity Vg.,.. 

^The sign of the product VgrJ determines the relative direc- 
tions between the current J and the group velocity Vgr- 



Note. — 



When n — regions III and IV disappear (see §4). 
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Table 2. Numerical solutions for gain and transmission for the shearing sheet model 



n 




K/flc 


Ivyh 






Amplification*^ 





1.0 


1.0 


0.3 


1.000496 


4.96E-4 


V 





1.0 


1.1 


0.4 


1.000447 


4.47E-4 


V 





0.8 


0.2 


0.3 


1.524428 


5.24E-1 


V 


1 


1.0 


1.0 


0.3 


3.89E-1 


4.60E-7 


X 


1 


1.0 


1.1 


0.4 


4.94E-1 


1.60E-6 


X 


1 


0.8 


0.2 


0.3 


9.38E-1 


7.63E-3 


X 



'^The gain, i.e., the refiection coefficient, defined by eq. (A13). 

^The transmission coefficient defined by eq. (A14). 

'^Amplification {G > 1) is marked with "y^' ; De-ampfification {G < 1) 
is marked with "x". 
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Fig. 1. — The potential U{r) = —kl{r) in the wave equation (28) for the case when n = 
0. There are two permitted regions and one forbidden region along the r-axis. The two 
permitted regions are: I (r < ri) and II (r > r2), where ri,2 are defined by a — ^k. The 
boundaries of the regions are marked by the two vertical dashed lines. The forbidden region 
is: A (ri < r < r2), which contains the corotation radius as marked by the dot. (When 
n = the corotation radius is not a singularity in the wave equation [28].) The directions of 
ingoing and outgoing waves are shown with horizontal arrows, which are defined relative to 
the corotation radius. The relations among the directions of the current, the wave-vector, 
and the group velocity in each permitted region are summarized in Table 1. 
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ingoing > 

outgoing < 



ingoing 
> outgoing 




Fig. 2. — The potential U{r) = —k'^{r) in the wave equation (28) for the case when n > 0. 
There is an intrinsic singularity at the corotation radius at r = Tc where a = and the 
potential is infinitely deep, as indicated by the cross-sign. There are four permitted regions 
and two forbidden regions along the r-axis. The permitted regions are: I (r < ri), II (r > r4), 
III (r2 < r < Tc), and IV (rc < r < r^), where ri^4 are defined by cr = ^ max {y/n ft ±, k), 
and r2,3 are defined by cr = ^ min {y/nfl±, k). The boundaries of the regions are marked by 
the four vertical dashed lines plus the f/-axis. The forbidden regions are: Al (ri < r < r2), 
and A2 {r^ < r < r^). The directions of ingoing and outgoing waves are shown with 
horizontal arrows, which arc defined relative to the corotation radius. The relations among 
the directions of the current, the wave- vector, and the group velocity in each permitted 
region are summarized in Table 1. 



